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AN EXPERIMENTAL AND ANALYTICAL VIBRATION STUDY OF 


ELLIPTICAL CYLINDRICAL SHELLS 

By John L. Sewall, William M. Thompson, Jr.,* 
and Christine G. Pusey 
Langley Research Center 

SUMMARY 

This paper reports an experimental and analytical vibration study of free-free 
isotropic elliptical cylindrical shells of constant mass ranging in cross-sectional eccen- 
tricity from zero (circular shell) to 0.916 (corresponding to a major-to-minor axis ratio 
of 2.5). Experimental resonant frequencies, nodal patterns, and mode shapes were 
obtained by use of an air-jet shaker or an electrodynamic shaker with a noncontact induc- 
tance probe that could be moved over most of the shell surface. Experimental frequencies 
are in generally good agreement with analytical frequencies calculated by means of a 
Rayleigh-Ritz type of analysis featuring multiterm longitudinal and circumferential modal 
expansions. As many as four longitudinal terms, including beam -vibration functions, 
were required along with an increasing number of circumferential trigonometric terms 
in each shell displacement series to obtain converged results with increasing eccentricity. 
As many as 32 circumferential terms were required for the highest eccentricity shell. 

Frequencies for shells with eccentricity of 0.916 were as much as 40 percent below 
the corresponding circular-shell frequencies. As eccentricity increased, fewer higher 
order experimental shell vibration modes could be identified, and there were more higher 
order analytical modes with two frequencies for each pair of longitudinal and circumfer- 
ential mode numbers. Analytical and experimental mode shapes for eccentricities of 
0.760 and 0.916 indicate considerable longitudinal and circumferential modal coupling that 
tends to preclude positive identification according to longitudinal and circumferential 
components. 

Frequencies calculated by methods involving various approximations to the noncir- 
cular curvature are in good agreement with both experimental and analytical frequencies 
of the present study for the inextensional modes, at least up to an eccentricity of 0.760 
(corresponding to a major-to-minor axis ratio of 1.538). 

Underwater Explosions Research Division, Naval Ship Research and Development 
Center, Portsmouth, Virginia. 



INTRODUCTION 


In the study of the response of shell structures to static and dynamic loading, con- 
siderably more attention has been given to shells of revolution than to shells of noncir- 
cular cross section. However, although many structural components of aerospace and 
submarine vehicles can be adequately treated as actual or equivalent shells of revolution, 
this approach may well be unacceptable for shells with significantly noncircular curva- 
ture. It is therefore profitable to examine the behavior of noncircular shells, not only to 
obtain basic data, but also to define the limitations of certain near-circular shell 
analyses. 

Some noncircular-shell investigations are reported in references 1 to 14. Refer- 
ences 1 to 5 are concerned with particular aspects of the static structural characteris- 
tics of slightly noncircular shells, and references 6 to 14 deal with dynamic characteris- 
tics. Specifically, references 1 to 3 are concerned with unstiffened and ring-stiffened 
oval cylinders subjected to hydrostatic pressure and show large effects of noncircularity 
on static stresses and displacements. Reference 4 presents a finite-difference analysis 
of stresses and displacements in an orthotropic noncircular shell of arbitrary curvature 
and varying circumferential thickness. Reference 5 presents series solutions to 
Donnell's equations for the displacements due to pressure loads for open noncircular 
shells with arbitrary end conditions. References 6 to 11 are vibration studies, of which 
references 7 to 9 indicate relatively small effects of noncircularity on the frequencies of 
free vibration of infinitely long oval cylinders. Reference 12 contains a considerable 
amount of experimental frequency and mode-shape data for clamped-free elliptical cyl- 
inders; however, all these cylinders had the same cross-sectional eccentricity. In ref- 
erence 14, the method of reference 5 is applied to the vibration of noncircular curved 
panels. Most of references 1 to 14 are entirely theoretical with no experimental data 
included for comparison with analytical results. Experimental and analytical static 
stresses are compared in reference 4. Reference 11 contains experimental frequencies 
for a variety of noncircular shells and includes comparisons of these frequencies with 
approximate analytical frequencies based on frequency equations of circular shells with 
equivalent radii. 

The purposes of the present paper are to compare measured and calculated vibra- 
tion data of short, free-free elliptical cylinders of constant mass and to determine the 
effects of cross-sectional eccentricity on shell frequencies and mode shapes. Attention 
is given to the effects of eccentricity (or noncircularity) both on the two sets of inexten- 
sional modes, as in references 9 and 10, and on two sets of higher order longitudinal 
modes. Vibration tests were conducted on four thin-shell isotropic cylinders of equal 
length, perimeter, and thickness, and with eccentricities ranging from zero (circular cyl- 
inder) to 0.916, corresponding to major-to-minor axis ratios from 1 to 2.5, respectively. 
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Analytical frequencies were obtained by application of an energy approach utilizing the 
Rayleigh-Ritz procedure in multiterm circumferential and longitudinal modal expansions. 
Trigonometric terms were used for the circumferential series and beam functions for 
the longitudinal series. Measured frequencies are compared with frequencies computed 
both by this analysis and, where possible, by the analyses of references 9 to 11. In ref- 
erences 9 and 10, the frequency is approximated by expressions involving the frequency 
of the circular cylinder plus power-series terms containing eccentricity parameters for 
ovals. The approximate frequency expressions of reference 11 are essentially circular- 
shell equations with equivalent or average values of the radius of curvature chosen to 
represent the noncircular cross section. 

Analytical mode shapes obtained in the present analysis are included to illustrate 
the nature of circumferential and longitudinal modal coupling due to noncircularity, and 
many of these mode shapes are compared with experimentally determined mode shapes 
for the radial (or normal) displacement. Pertinent details of the analysis are given in 
appendix A. In addition, analytical frequencies of the elliptical shell models with 
assumed freely supported ends are included in appendix B, and comparisons are made 
with the analytical frequencies of reference 10. 

SYMBOLS 


^ijmp»®ijmp^ijmp | 

■ ijmp’GijmpjHjjmp 


elements of stiffness matrix (see appendix A) 


a,b length of semimajor and semiminor axis, respectively, of elliptical shell 


a mp>kmp> c mp generalized coordinate for antisymmetric-mode displacement u, v, 

and w, respectively 

%p>^mp>%p generalized coordinate for symmetric-mode displacement u, v, 

and w, respectively 


E Young's modulus 


e 


eccentricity of elliptical shell, 



f circular frequency, hertz 


h shell thickness 
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Il to Ig 
i 

J 1 to J 16 

j,P 

K 

L 

M 

m 

N 

n 

P,Q 


v P < l > 

R 

R 0 

s 

s o 

T 

t 

U 


circumferential integrals (see appendix A) 
longitudinal mode number; \J^1 in equations (5) 
longitudinal integrals (see appendix A) 

half the number of longitudinal node lines for elliptical shells 
stiffness submatrix (see eq. (9)) 
length of shell 

mass submatrix (see eq. (9)) 
longitudinal mode number 
beam eigenvalue (see appendix A) 
circumferential mode number for circular shell 

integers identifying upper limits on number of modes used in the analysis 
(see eqs. (5)) 

generalized coordinate 

radius of curvature of shell 

arbitrary reference radius 

arc length around shell 

total circumference of shell 

kinetic energy 

time 

strain energy 
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u,v,w 


longitudinal, circumferential, and radial shell displacement, respectively 
(see fig. 6) 

x mu> x mv> x mw mth longitudinal mode shape for u-, v-, and w-displacement, respectively 
x,y,z longitudinal and cross-sectional coordinates (see fig. 6) 

“ijmpj&jmpj'l'ijmp elements of mass matrix (see appendix A) 
y beam eigenvalue property (see appendix A) 

6x, £$,6x0 membrane strains (see eqs. (1)) 

6 angle between the minor axis and the radius vector R (see fig. 6) 

K x’ K 0’^x0 changes of curvature (see eqs. (1)) 
p. Poisson’s ratio 

p shell mass density 

<i> = 1 - e^cos^p 

<p central angle (see fig. 6) 

c o angular frequency, radians/second 

Subscripts: 

a antisymmetric 

i, m identifies ith and mth longitudinal modal components 

j, p identifies jth and pth circumferential modal components 

s symmetric 

Superscript T denotes transpose of a matrix. 

Matrix elements and integrals with bars are associated with symmetric mode. 

5 

I 



I Minim ■■mill ■■■ nan hi iih 


■ hi in Minn i ■■■ ■mm 


Dots over quantities denote differentiation with respect to time. 
Primes denote differentiation with respect to x, for example, = 


®mW 

dx 


d 2 X 


mw 


(x) 


k mw 


cbr 


and 


EXPERIMENTAL INVESTIGATION 
Models 

Four cylindrical shell models of equal perimeter, length, and thickness were used 
in the investigation. One model was circular in cross section, and three were elliptical. 
Each model was 24 inches (61.0 cm) long and was made from two shaped sections butt- 
welded together with smoothed seams at the ends of the major axis. Each section was 
6061 aluminum alloy and was 0.032 inch (0.813 mm) thick. Keeping the perimeter, 
length, and thickness the same for all models maintained a constant total mass. The 
model dimensions and eccentricities, along with the material properties E, p, and p, 
are indicated in table I. 


Test Apparatus and Procedure 

Vibration tests were conducted with each model suspended by 12 soft elastic sup- 
ports attached at equidistant points around the perimeter at one end of the cylinder, as 
shown in figure 1. The basic parts of the test apparatus are the same as those reported 
in reference 15 and are shown schematically in figure 2 as arranged for some of the tests 
in this investigation. Each model was excited by an oscillatory force directed normal to 
the shell surface. Two types of vibration exciters were used, an air-jet shaker of the 
type described in reference 16 and a small electrodynamic shaker capable of a 1-i- pound 

(6.67-N) vector force output. At low frequencies (less than 50 Hz), the air shaker pro- 
vided a noncontact excitation force which permitted unrestricted, large-amplitude motion 
of the cylinder wall. For relatively low-amplitude responses at higher frequencies, the 
electrodynamic shaker was used and was attached to the cylinder wall by means of a 
lightweight vacuum cup. Of the two shakers, the air shaker was believed to give more 
reliable lower frequencies because of the appreciable stiffness and mass effects of the 
electrodynamic shaker on these frequencies. (See also ref. 17.) Stationary and movable 
noncontact inductance probes of the type described in reference 15 were used to determine 
resonant peaks, phase shifts, and mode shapes in normal displacement. Resonant fre- 
quencies were obtained by tuning the shaker for maximum amplitude response normal to 
the shell wall at an antinode. 
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With the test setup of figure 2, nodes of the normal displacement were located by 
monitoring the outputs of the fixed and movable probes in the forms of Lissajous figures 
on an oscilloscope and by observing the phase shifts for different positions of the movable 
probe. No mode-shape measurements were made. This apparatus was used for all the 
tests on model 2 and for some of the tests on model 4. 

For the rest of the test program, normal mode shapes were measured with the aid 
of the motorized movable probe apparatus shown in figures 1 and 3. The movable induc- 
tance probe was driven by an electric motor along a track mounted on a concentric, 
3/8-inch (0.952-cm) steel elliptical strap surrounding the model, as shown in figure 1, 
in order to obtain circumferential mode-shape components. This strap was supported on 
four vertical threaded rods and could be rotated by another electric motor through a 
chain drive in order to obtain longitudinal mode-shape components. The supporting 
instrumentation was more involved than that represented in figure 2 and is detailed in 
reference 15. Figure 3 is a closeup view of the motorized probe assembly between the 
elliptical track support and the model. The thin white wire leading upward from the back 
of the probe itself goes to the electronic instrumentation which controls the distance of 
the probe from the model and the operation of an x-y plotter that traces out normal shell 
displacements as the probe travels over the shell surface. 

Experimental Results 

Resonant frequencies of the four cylindrical shells are listed in tables II to IV. 

The mode identification integers m, n, and p used in tables II and in are defined in 
figure 4. Mode shapes corresponding to frequencies in tables in(b) and IV for model 3 

(|= = 1.538, e = 0.76oj are shown in figure 5. The air shaker was used for frequencies 

up to about 170 Hz for the circular shell and up to about 50 Hz for the elliptical shells, 
and the electrodynamic shaker (designated as "electric shaker" in the tables and figures) 
was used for all other modes (except as noted in tables II and ni(c)). Erratic random 
deviations of some of the mode shapes in figure 5 were due to occasional drifts from 
resonance during mode-shape measurement. This drift was more apt to occur with the 
air shaker than with the electrodynamic shaker. The closely spaced ripples superim- 
posed on the responses at the lowest frequencies are attributed to the closeness of these 
frequencies to the natural frequency of the probe servo system. 

The modal identifications in figure 4 apply to normal (or radial) shell displace- 
ments. The longitudinal mode number m denotes the number of circumferential nodes. 
The circumferential mode numbers n and p denote, respectively, the number of cir- 
cumferential waves for the circular shell and one-half the number of longitudinal nodes 
for the elliptical shells. The modes for m = 0 and 1 are designated as "inextensional" 
because of their straight-line variations in the longitudinal direction and because of the 
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monotonically increasing frequencies with increasing n or p, consistent with purely 
bending vibrations of the shell. Modes for m = 2 are characterized by curved (indi- 
cating elastic) longitudinal mode-shape components and by the existence of relative mini- 
mums in the variation of frequency with n or p consistent with the well-known shell 
behavior associated with combined bending and stretching (or membrane) deformations. 

Nearly all resonances excited for the circular shell (table n) and most of the lower 
frequency inextensional modes of the elliptical shells (table Ht) could be clearly classi- 
fied according to figure 4. In addition, the identifiable elliptical-shell modes were either 
symmetric or antisymmetric with respect to the minor axis, and this distinction is 
clearly evident in figure 5 for the lower frequency modes. The symmetric modes were 
excited by alining the shaking force along the minor axis, and the antisymmetric modes 
were obtained with the shaker located part way between the ends of the major and minor 
axes at a point of relatively large displacement. Most symmetric and antisymmetric 
frequencies differed bj less than 1 Hz. 

The number of identifiable resonances decreased with increasing frequency and 
increasing eccentricity e until few elliptical -shell modes could be exclusively identi- 
fied as m = 2or 3. Many of these modes were highly coupled with m = Oor 1 and/or 
more than one p mode, as is evident in figure 5 by the different values of m and/or 
p in different parts of the model. Moreover, both symmetric and antisymmetric char- 
acteristics were present in the circumferential mode shapes of most of these modes. 
Localized effects due to the electrodynamic shaker also contributed to the uncertainty 
and ambiguity of modal identification, mainly by causing local longitudinal mode-shape 
distortions that tended to become more severe with increasing frequency, as is shown in 
figure 5. Consequently, the frequencies in table IV are simply listed in order of 
increasing value. 


METHOD OF ANALYSIS 

The analytical method derived and applied in this paper is an application of the 
well-known Rayleigh-Ritz procedure. The assumed forms for the inplane and radial (or 
normal) displacements were finite series representing circumferential and longitudinal 
(or axial) components of these displacements. Elementary trigonometric functions rep- 
resented the circumferential components for both symmetric- and antisymmetric-mode 
frequencies. Longitudinal displacement components were chosen to satisfy geometric 
end conditions (displacement and slope). For the free-free end conditions, longitudinal 
mode-shape components for the two lowest sets of modes (m = 0 and 1) were approxi- 
mated by simple constant and linear algebraic expressions, and as in reference 18, the 
longitudinal mode-shape components for the higher modes (m = 2) were approximated by 
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free-free beam-vibration functions. Longitudinal mode-shape components are also 
included for the freely supported end conditions and are given by trigonometric functions. 


Strain-Displacement Relations 


The elliptical geometry embodied in the analysis is detailed in figure 6, and the 
analysis is formulated on the basis of thin-shell theory in terms of circumferential and 
axial coordinates 0 and x. The deformation of the shell middle surface is expressed 
in terms of the displacements u, v, and w in the axial, circumferential, and normal 
directions, respectively, with the positive normal direction considered to be outward. 

By use of Sanders' theory (ref. 19), the strains eg, and e x g and changes of curva- 

ture k x , Kg, and H^g are given in terms of the displacements by the relations 




9u 


£x ' 9x 


1 /9v 

e o-R[Je + w 


, _ l/Bv 1 9u 

2\9x R 90, 


K x =~ 


9 2 W 

9x 2 


1 

9 2 W 

9v 

1 

dR 

h - v\ 

R 2 

90 2 

90 

R 

d0 ( 


9v _ 1 

. a 2 w 


1 

9u 



( 1 ) 


xd 4R 9x " R 9x 90 " 4r 2 90 




where R is the radius of curvature which, for the elliptical cross section, is given by 


R = 


a $ 3 / 2 


( 2 ) 


with = 1 - e 2 cos 2 0 and 
larity of the cross section. 



The 


dR 

d0 


term of 


Kg is due to the noncircu- 


Strain and Kinetic Energies 

The strain energy of an isotropic shell of uniform thickness h may be written as 
follows in terms of the strains and changes of curvature: 


U = 


Eh 


2(l - fi) 


J 0 Jq 4 + 2 ^x e 0 + e f + 2 ( 1 - ^ 0 ; 


R d0 dx 


f^2) C C [** + 2lM<xKo + k2 o + 2 ( 1 - 


24(: 

where E is Young's modulus and p is Poisson's ratio. 


R d0 dx 


( 3 ) 
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The shell kinetic energy is 


T 


ph 

2 



+ v c 


+ 



R d6> dx 


where p is the mass density of the shell. 


(4) 


Modal Functions 


Each displacement u, v, or w is assumed to be a finite series, each term of 
which is the product of circumferential and axial modal functions weighted by an appro- 
priate amplitude coefficient (or generalized coordinate) q mp (t). With the assumption of 
simple harmonic motion of frequency u>, the assumed displacements may be written as 

P Q 

u(x,0,t) = J ^ (a mp cos p0 + a mp sin p6>)x mu (x)e lwt 
m=0 p=0 


p Q 

v(x,0,t) = J 2 (^mp sin PS - b mp cos P#)x mv (x)e lwt \ 
m=0 p=0 j 


(5) 


P Q 

w(x,e,t) = ^ y (c mp cos P0 + c mp sin pejx mw (x)e lwt 
m=0 p=0 


where a mp , b mp , and c mp are amplitude coefficients for the symmetric modes and 
a mp , b mp , and c mp are amplitude coefficients for the antisymmetric modes. Upper 
limits for the series of axial and circumferential modal functions are denoted by P and 
Q, respectively. 

The X m functions in equations (5) approximate the longitudinal modal components 
which are chosen to satisfy displacement and slope conditions at the ends of the shell. 

For example, the conditions of zero displacement for v and w at the ends of a simply 
(or freely) supported shell without axial constraint are satisfied by the familiar trigono- 
metric functions 


x mu 

x mv 


= cos 


m7?x 

L 


= x mw = sin 




imrx 


L 


J 


( 6 ) 
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The arbitrary end displacements and slopes of a free-free shell can be satisfied by 
choosing the following modal functions: 


x 0u = o 

V . V 

A 0v “ J *"0w 



X 



■^lv “ ^lw 



(m = 0) (7a) 


(m = 1) (7b) 


^mu 

V 

A mv 


= xl 


m 


X mw = X m = cosh N m x + cos N m x )> 
- y m ( sinh N m x + sin N m xj 


(m S 2) (7c) 


where X m = — — — and N m and y m are the eigenvalue properties of a free-free 

vibrating beam, as tabulated in reference 20, for example. The functions in equa- 
tions (7c) are essentially the same as those used in reference 18. The existence of the 
inextensional mode shapes approximated by equations (7a) and (7b) is demonstrated 
experimentally in this paper and, for free-free circular cylinders, in reference 21. 


Derivation of Frequency Equation 

With the substitution of equations (1), (5), and (7) into equations (3) and (4), the 
equations of motion are obtained from the following relations consistent with the 
Rayleigh-Ritz procedure: 


-!- u(x,e) - u 2 T (x,e)l =-|— [u(x,e) - u> 2 T (x,e)] = -*- 
da ijL j ab ij L J BCij 


T— — |U(x,0) - w 2 T(x,e) 
a ij L 


9bij 


u(x,e) - co 2 t(x 


, 9)1 = 


8C ii 


|U(x,0) - w 2 T(x,0) 


u(x,0) - u> 2 t ( x,e) 


= o 


= 0 


( 8 ) 
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The operations represented by equations (8) lead to the familiar eigenvalue-eigenvector 
formulation which may be expressed in the following general matrix form: 



where K and M represent square stiffness and mass matrices of size 3(P + 1) 
by Q + 1, q denotes an eigenvector, and the subscripts s and a identify sym- 
metric and antisymmetric matrix elements, respectively. The eigenvalue is given by 

u ) 2 r 2 p (1 - fj2) 

A = . 

E 

The off-diagonal blocks vanish, and equation (9) uncouples into two simpler equa- 
tions, one for the symmetric and the other for the antisymmetric modes. The equation 
for the symmetric mode may be put in the form 



The equation for the antisymmetric mode is identical except the matrix elements do not 
have bars. The superscript T denotes the transpose of a matrix. Each letter in the 
square matrices of equation (10) represents a submatrix, the elements of which are given 
in detail in appendix A. Longitudinal modal components of the matrix elements for the 
freely supported shell, based on equations (6), are given in appendix B. 

For the circular cylindrical shell (e = 0), R is constant and can be considered to 
be equal to R 0 , and each series of circumferential terms in equations (5) reduces to a 
single trigonometric function. Equation (10) then gives the same frequencies for sym- 
metric and antisymmetric modes. 

ANALYTICAL RESULTS AND COMPARISON WITH EXPERIMENT 

The analysis described was applied to the circular and elliptical cylindrical shells 
of table I. Analytical frequencies of the free-free shells are included along with experi- 
mental frequencies in tables II to IV. Analytical mode shapes are presented in figures 7 
to 12, and experimental mode shapes are included for comparison in figures 7 and 9 to 11. 
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Circumferential mode-shape components are presented in figures 7 to 10 and 12, and 
experimental components included in figure 9 are the topmost mode shapes shown for 
each frequency in figure 5. Longitudinal mode-shape components are presented in fig- 
ure 11 for model 4. Experimental nodes obtained with the simplified test apparatus of 
figure 2 are included in figure 10(a) for model 4. Other frequency comparisons are made 
in figures 13 and 14. 

Elliptical-shell frequencies calculated by the method of the present paper are also 
compared in table III with frequencies calculated by the method of references 9 to 11. 
Analytical frequencies of the freely supported shells calculated by the present analysis 
are given in appendix B and are compared with analytical frequencies of reference 10 in 
figures 15 and 16. 


Convergence of Analytical Results 

Shell modes . - In the following table are listed the numbers of the longitudinal and 
circumferential functions used in each of the displacement series of equations (5) to 
obtain converged analytical frequencies and mode shapes: 


Model 

Tables 

e 

Number of terms in each series of 
equations (5) 




Longitudinal 

Circumferential 

1 

II 

0 

4 

a l 

2 

IH(a), IV 

.526 

4 

13 

3 

m(b), IV 

.760 

4 

16 

4 

m(c), iv 

.916 

2 

32 


a p = n. 


The numbers in this table apply to all even-numbered or all odd-numbered modal func- 
tions, since, as noted in appendix A, even-numbered and odd-numbered functions do not 
couple when included together in equation (10). Results were considered to be converged 
when the analytical frequencies of tables II to IV were negligibly changed with the inclu- 
sion of more terms in the series than those in the foregoing table. Beginning with the 
circular shell (model 1), which required only a single circumferential term in each 
series, converged solutions for this model were obtained with four terms, including 
beam -vibration functions, in each longitudinal displacement series, that is, m = 0, 2, 4, 6 
and m = 1, 3, 5, 7 in separate calculations for each circumferential wave number n. 
(The circular-shell results in table II and in fig. 13 are actually based on five longitudinal 
modal functions, the fifth term being the beam functions for m = 8 and 9 in the foregoing 
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m sequences.) Analytical circular-shell frequencies of the inextensional modes 
(m = 0 and 1) were negligibly different from those in table II when just a single longi- 
tudinal modal function (m = 0 or 1) was included in each displacement series. 

Next, for the elliptical shells in order of increasing eccentricity, the number of 
longitudinal terms was held to four, and the number of circumferential terms was 
increased to achieve convergence. As the foregoing table shows, more circumferential 
terms were required with increasing eccentricity, and at e = 0.916 the number of 
required terms exceeded the capacity of the computing facility utilized in this study 
(Control Data 6600 computer system at the Langley Research Center). Hence, the num- 
ber of longitudinal terms had to be reduced to allow for more circumferential terms, and 
the combination of two longitudinal and 32 circumferential functions listed in the table 
gave results that are as close to convergence as could be obtained for this model. This 
observation is further borne out in comparisons of the analytical frequencies in 
tables HI(c) and IV with those in table V for model 4. Comparison of tables IH(c) with V 
shows frequencies for m = 0 and 1 to be satisfactorily converged with a single longi- 
tudinal term and 30 circumferential terms. However, it appears from tables IV and V 
that additional longitudinal terms were needed more than additional circumferential 
terms to obtain converged results for m = 2 and 3; and 32 circumferential terms appear 
to be enough. 

Analytical circumferential mode shapes are shown for various modal approxima- 
tions in figure 12. These mode shapes, together with their corresponding frequencies in 
table V, show the sensitivity of certain modes to variations in the number of circumfer- 
ential terms included in each displacement series of equations (5). The existence of dual 
frequencies at some values of p and none at others also occurred for models 2 and 3, 
as may be seen in figures 8 and 9(e) to 9(h), and there appears to be no consistency to 
this behavior other than the increasing number of such modes with increasing eccentricity. 

Rigid-body modes . - The lowest frequencies of free-free shells are zero frequencies 
for six rigid-body modes. These zero frequencies occur at p (or n) = 0 and 1 for the 
inextensional modal families m = 0 and 1, and in table VI are listed the displacements 
with zero frequencies for both circular and elliptical shells. Six zeros were obtained for 
the circular shell regardless of the number of terms in the displacement series, and six 
zeros were obtained for two of the elliptical shells (e = 0.526 and 0.760) with the same 
number of terms as needed for convergence of the shell modes. As shown in the last 
column of the table, five zero frequencies were obtained for the remaining high- 
eccentricity elliptical shell (e = 0.916), but the sixth frequency was very low, 1.01 Hz. 

With only a single modal product in each displacement series, only three rigid-body 
modes could be obtained for any of the elliptical shells. 


14 



Effects of Eccentricity 


The close agreement between analytical and experimental frequencies and mode 
shapes of the circular cylindrical shell (model 1) in table II and figure 7 afforded a satis- 
factory basis for evaluating eccentricity effects of the elliptical shells. Both analytical 
and experimental frequencies and the corresponding circumferential mode shapes of the 
elliptical shells showed greater distortions from sinusoidal-type responses as eccentric- 
ity increased. Tables in and IV show generally good frequency agreement over the 
eccentricity range of the models, and corresponding analytical and experimental mode 
shapes in figures 9 and 10 are fairly well correlated. This correlation is clearly better 
for the low-frequency inextensional modes than it is for the higher frequency modes, 
particularly those involving m = 2 and 3. Experimental mode shapes are compared, 
somewhat arbitrarily, with analytical symmetric or antisymmetric mode shapes 
according to shaker position as indicated in figure 9, although, as previously noted, most 
modes involving m = 2 and 3 had both symmetric and antisymmetric characteristics 
regardless of shaker position. For the experimental mode shapes included in fig- 
ure 10(b), the shaker was located at the end of the minor axis. 

Variations of frequency with mode numbers are shown in figure 13 for eccentric- 
ities of zero (model 1) and 0.526 (model 2). Similar plots for the higher eccentricity 
shells were found to be impracticable because of the dual frequency behavior noted 
earlier. This property for model 2 is identified in the right-hand plot of figure 13 by 
the missing analytical frequency for m = 3, p = 10 and the second analytical frequency 
at m = 3, p = 12 (denoted byB). 

Eccentricity effects on both experimental and analytical frequencies are summa- 
rized in figure 14. Except for the experimental inextensional frequencies for e = 0.526 
(table 111(a)), a general decrease in frequency occurred with increasing eccentricity, and 
this decrease was largest for e § 0.760 (or a/b = 1.538). This effect was much 
smaller for the inextensional modes than for the minimum frequencies of m = 2 and 3 
which, at e = 0.916, were approximately 40 percent lower than the minimum circular- 
shell frequencies. For e = 0.526, there was a slight, though consistent, rise in most 
experimental inextensional frequencies of about 2 or 3 percent, which was not predicted 
by the analysis. 


Comparisons With Other Methods 

Inextensional frequencies calculated by the methods of references 9 to 11 are 
included in table III and appear to be in almost as good agreement with experimental fre- 
quencies as are the frequencies calculated by the present analysis, at least for e = 0.760. 
The methods of references 9 and 10 are based on approximations to the noncircular 
radius of curvature that are contained in the terms of power series which include the 
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circular-shell frequency as the first term. The method of reference 9 is valid for 
a/b < 2.06 and is thus not applicable to e = 0.916 (a/b = 2.5) in the present study. 

The noncircular terms of reference 10 are based on a Fourier series representation of 
the radius of curvature in which the coefficients of the series may be obtained by approx- 
imating the noncircular cross section by two pairs of circular arcs, each having a differ- 
ent radius. By following the recommendations of both of these references, the frequen- 
cies in the two right columns of table in were obtained by retaining only a single 
noncircular power-series term, other than the circular-shell frequency, in the frequency 
equation. Most of the frequencies calculated in reference 9 are seen to be somewhat 
higher than those calculated in reference 10, the differences for e = 0.760 being greater 
than those for e = 0.526. (Compare tables HI(a) and (b).) 

The applicable frequency equation of reference 11 is essentially the circular-shell 
equation for m = 1 with the radius of curvature approximated by an average value 
based on the semimajor and semiminor axes of the elliptical shell. The frequencies 
given by this method are in about as good agreement with experiment as those of refer- 
ences 9 and 10 for e = 0.526 and 0.760 and tend to be slightly higher than those of the 
present analysis and reference 10 at higher values of p. However, the agreement at 
e = 0.916 is clearly poorest of all the methods considered, the discrepancy being as 
much as 14 percent (at p = 5). 


C ONC LUDING RE MARKS 

This paper reports a vibration study of free-free aluminum elliptical cylindrical 
shells of constant mass over a range of cross-sectional eccentricities from zero (cir- 
cular shell) to 0.916 (corresponding to a major-to-minor axis ratio of 2.5). With 
increased eccentricity over this range, shell frequencies were reduced by as much as 
40 percent. Experimental resonant frequencies, normal-displacement node locations, 
and mode shapes were obtained by use of an air-jet shaker or an electrodynamic shaker 
with a noncontact inductance proximity sensor that could be moved over most of the shell 
surface. Experimental frequencies are in generally good agreement with analytical fre- 
quencies calculated by means of a Rayleigh-Ritz type of analysis featuring multiterm 
longitudinal and circumferential modal expansions. Single longitudinal functions in each 
of the shell displacement series were sufficient to obtain converged results for the two 
families of inextensional free-free modes (longitudinal mode numbers of 0 and 1). Four 
longitudinal terms, including beam -vibration functions, were required in the series to 
obtain converged results for the next two higher order longitudinal modes (longitudinal 
mode numbers of 2 and 3). Along with these longitudinal approximations, an increasing 
number of circumferential trigonometric terms were required for convergence with 
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increasing eccentricity, and as many as 32 circumferential terms were required for the 
highest eccentricity shell. 

As eccentricity increased, fewer higher order experimental shell vibration modes 
could be identified exclusively according to longitudinal and circumferential mode num- 
bers, and there were more higher order analytical modes having two frequencies for 
each pair of longitudinal and circumferential mode numbers. Analytical circumferential 
mode shapes shed further light on this modal coupling and dual frequency behavior by 
indicating the presence of more than one longitudinal and/or circumferential mode in 
each mode shape. 

Frequencies calculated by methods involving various approximations to the non- 
circular curvature are in good agreement with both experimental and analytical inexten- 
sional frequencies of the present study for eccentricities up to 0.760 (corresponding to a 
major -to -minor axis ratio of 1.538). 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va. , September 30, 1970. 
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APPENDIX A 

MATRIX ELEMENTS AND INTEGRALS IN RAYLEIGH-RITZ 
VIBRATION ANALYSIS 

This appendix contains detailed expressions of the matrix elements in equation (10). 
Also included are the circumferential and longitudinal integrals involved in these ele- 
ments. The circumferential integrals were evaluated numerically by use of a 10-point 
Gaussian quadrature method with 20 subintervals for a total of 200 integrating stations. 

The matrix elements in equation (10) may be written as follows for the symmetric 
modes of an elliptical cylindrical shell: 


2 _ 


Aijmp = f Vl + ^ IP ' f 2 J 2 + fefiP f 3 J 2 


(Ala) 


®ijmp ~ *6 jp P|iJ3 " 1 


; i - a 


J 4 + — n 

2 lb 2 . 


a \ 2 jh 2 (l - n) f 


32 


I 4 J 4 


(Alb) 


2 2 

C ijmp - ^6jp(l + 60j50p)j 5 + 3P h U 24 ~ M) I 4 J 6 


(Ale) 


Fiimo = ^ 3P I 5 J 7 + Hr ^ J 6 J 8 + l~) ^ + ^ * 2 J 8 


ijmp b2 


b 2 1 - p. 


'a \ 3 h 2 


a 3h 2 (l - p) 7 . 


a 2 


b 


(Aid) 


a -f T a h' 


a\ h 2 n f„ r n x a ph 2 (_l - M) T ’ 


^mp^iVo-^^lO + J 12 ^ ^9 


(Ale) 




_a h 2 
b 2 


_a \ h 2 
Vb 2 




2 d 16 


(A If) 
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APPENDIX A — Continued 


OL • • 

ijmp 


l“l J 2 

sRi 


(A2a) 


^ijmp " 



(A2b) 


r ijmp 



(A2c) 


where I is a circumferential integral and J is a longitudinal integral and where 
6jp, 5 0j , and 8 q p denote Kronecker delta functions of the form 


61d '{l 

(k = l) 
(k * 1) 


k,l = 0, 1, 2, . . j,p. 



The circumferential integrals are defined as follows: 



Ij = T 3>~3 /^cos j© cos p© d© 
J 0 


(A3 a) 

12 = T $^/^sin j© sin p© d 0 
J 0 


(A3b) 

Io = C t&^/^sin j© sin p© d© 
J 0 


(A3c) 

n2n o 

1a = \ $ sin j© sin p© d© 

J 0 


(A3d) 

Ic = C <3?^/^cos j© cos p© d© 
J 0 


(A3e) 

Ig = C f,"3/2sin jg S j n pg ,jg 


(A3f) 


0 
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APPENDIX A — Continued 


I 7 = ^ cos ]6 + e^$"^sin 26 sin jdj|p cos P 0 + e 2 $"*sin 26 sin pfljdfl 


fg = ^ cos j£ + ^ e ;s $ - ' l sin 26 sin jfijcos p 0 dd 


2 q2a" ^ci 1 


fg = y <3>“^cos j 0 ^p cos p6 + e 2 $“^sin 2 6 sin j^jd# 


(A3g) 

(A3h) 

(A3i) 


The antisymmetric matrix elements are identical in form to equations (Al) and (A2) and 
are obtained by simply interchanging sin j0, sin p6 for cos j 0, cos p 0 inequa- 
tions (A3). In addition, 

2 2 

Cijmp = *M6jp(l - 5 0 jS0p]j5 + (^j iP h { \~ M) Ve (A4) 

For both symmetric and antisymmetric matrix elements, it can be shown that the integrals 
in equations (A3) vanish if j ± p is odd and have values if j ± p is even (i.e., 
j,p = 0, 2, 4, . . .or 1, 3, 5, . . .). Thus, two calculations must be made in every appli- 
cation of the analysis for symmetric or antisymmetric modes, one for all even circum- 
ferential modal functions and another for all odd functions. 


Longitudinal Integrals for Arbitrary End Conditions 


The longitudinal integrals are defined by the following general expressions, which 
can be evaluated in closed form for particular choices of modal functions and end 
conditions: 


J 1 4 0 L X iu X mu 

r L 

J 2 - \ x iu x mu dx 
J 0 


J 5 

J 6 


io 

r L 

x iu x; 

n 


mw 


mw 


dx 


dx I 


(A5) 


J3 

J 4 



x iu x mv 


X iu x mv 


dx 


dx 



(Equations continued on next page) 
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APPENDIX A - Continued 


*-r 

J o 
J o 


x iv x mw dx 


d 10 X iv x mw ^ 


J 11 = x !v x mw dx 


'12 


- T x iw x mw dx 


d 13 " X iw X mw dx 

r L 

J 14 = J Q X iw X mw dx 

r L 

d 15 “ jp 

- r 

^0 


"mw 


(A 5) 


Longitudinal Integrals for Free-Free Shell 

With the ends of the cylindrical shell free and the longitudinal modal functions given 
by equations (7), the integrals in equations (A5) reduce to the following forms on the basis 
of reference 22: 

i = 0 (eq. (7a)): 

J — J ^ = J ^ = = J ^ — J 0 = J g = J ~ J = *^2.5 = *^10 = ^ (All v3.1u.gs of m) 

0 


J 7 - Jg - J i2 ~ 


J 10 = J 14 = 


2N m y m [l - (-l) m+1 ] 


0 


(i * m) 

(i = m)|> (A6a) 

(m g 2) 

(m = 0,1) 

J 


i - 1 (eq. (7b)): 

J 1 = J3=J 5 = J 1 3 = J 15 = 0 


J 2 - J 4 - J 6 - J 8 _ J 11 ~ J 16 _ / 


r 

0 

1_ 

L 


2_ 

L 


1 + (- 1 ) 


m+1 


J 


(All values of m) 

(m = 0) 

(i = m) 

(m S 2) 


) (A6b) 


(Equations continued on next page) 
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APPENDIX A — Concluded 


Jrj - Jg - J 12 - \ L 
12 


2Nm7m 


J 10 _ J 14 _ ( 0 


1 _ (_i)m+l 




i,m S 2 (eq. (7c)): 


J 1 = J 13 


N^L 
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1N m 1N i 


N mr m y m N m L + 6 
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4N i ^m N m " n N i 


J 3 = J 5 = J 15 = \ N m “ N i 

Nm7m(2 - 7m^ m^) 


(-l) 1+m + 1 




J 7 = J 9 = J 12 = 


4Nm(nNj - 7m% 


J 10 _ J 14 “ \ 
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Ni "Nm 


N m"/m( 2 “ 7 m N m L 


(- l) 1+m + 1 


(i * m) 
(i = m) 

(m 12) 
(i = m) 
(m = 0) 


(i * m) 
(i = m) 

(i * m) 

(i = m) 

(i * m) 

(i = m) 

(i * m) 
(i = m) 

(i * m) 

(i = m) 

J 


> (A6b) 


> (A6c) 


The coupling characteristics noted in the body of the paper are implicit in equa- 
tions (A6), namely, the integrals vanish if i ± m is odd and have values if i ± m is 
even (i.e., i,m = 0, 2, 4, . . .or 1, 3,5,. . .). 
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APPENDIX B 


APPLICATION OF ANALYSIS TO FREELY SUPPORTED 
ELLIPTICAL SHELLS 

This appendix demonstrates the application of the elliptical-shell vibration analysis 
developed in the body of the paper to the models of table I with assumed freely supported 
end conditions. With the use of equations (6), each series of longitudinal terms in equa- 
tions (5) reduces to a single trigonometric function. The integrals in equations (A5) used 
in equations (Al) and (A2) thus reduce to the following for this end condition: 


Jl = J 8 = Jii = Ji6=^— 

(Bla) 

j 2 = j 7 = j 9 = j 12 

(Bib) 

T T m77 

J3 = J 5 = ~j- 

(Blc) 

T , m,r 

J 4 = J 6 = — 

(Bid) 

T T T ( m ^) 2 

J 10 = J 14- J 15 " 2L 

(Ble) 

Jl3 = (m f 

OT *5 

(Blf) 


2L" 3 


where, for this end condition, m is defined as the number of longitudinal (or axial) half- 
waves. The reduction of each series of longitudinal terms in equations (5) to a single 
trigonometric function is due to the relations 

C L sin H* sin dx = f L cos cos ™ dx = 0 (i * m) 

dq L L dq L L 

By use of equations (Bl) and equations (6) in equation (10), calculated frequencies 
were obtained for the models of table I with the same numbers of circumferential terms 
as those used for the free-free models. Results are given in tables VII and VIII and are 
compared in figures 15 and 16 with frequencies calculated by reference 10. The two sets 
of results are seen to be in close agreement, and the decrease in frequency with increasing 
eccentricity is similar to that found for the free-free shells for m = 2 and 3. 
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TABLE I.- PROPERTIES OF CYLINDRICAL SHELLS 

|e = 10 7 lb/in2 (68.95 GN/m2); M = 0.3; p = 2.588 x 10~ 4 lb-sec2/in 4 (2768 kg/m3); 
L = 24 inches (61.0 cm); h = 0.032 inch (0.813 mm)J 



. -- 

Model 

a 

b 

a/b 

e 

1 

12.00 in. 
(30.50 cm) 

12.00 in. 
(30. 50 cm) 

1.0 

0 

2 

12.95 in. 
(32.90 cm) 

11.01 in. 
(27.97 cm) 

1.176 

0.526 

3 

14.39 in. 
(36.55 cm) 

9.35 in. 
(23.75 cm) 

1.538 

0.760 

4 

16.39 in. 
(41.63 cm) 

6. 56 in. 
(16.66 cm) 

2.50 

0.916 
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TABLE n.~ EXPERIMENTAL AND ANALYTICAL FREQUENCIES 
OF FREE-FREE CIRCULAR CYLINDRICAL SHELL 






Frequency, Hz 




n 

n 

i = 0 

m = 1 

m = 2 

m = 3 


Analysis 

Experiment 

Analysis 

Experiment 

Analysis 

Experiment 

Analysis 

Experiment 

0 





a 2733.0 


a 2594.0 


1 

0 


a 2014.00 


a 2293.0 


a 2494.0 


2 

5.645 

5.6 

7.54 

7.7 

1616.0 


2101.0 


3 

16.0 

15.6 

19.0 

18.9 

1068.0 


1709.0 


4 

30.6 

30.0 

34.2 

35.7 

717.8 


1340.0 


5 

49.5 

48.9 

53.4 

53.0 

504.8 


1045.0 


6 

72.6 

72.0 

76.7 

76.4 

375.6 

b 377.3 

823.5 


7 

100.0 

99.3 

104.1 

103.8 

299.9 

299.1 

663.4 


8 

131.5 

131.0 

135.7 

135.3 

262.2 

c 257.4 

551.4 








262.1 



9 

167.25 

166.9 

171.5 

d 170.7 

253.4 

e 248.8 

477.8 

484.4 



f 167.2 




249.3 



10 

207.2 

206.9 

211.5 

210.2 

266.3 

268.8 

436.3 


11 

251.4 

250.5 

255.7 

253.0 

294.7 

290.9 

421.8 

422.1 









425.6 

12 

299.7 

e 301.6 

304.1 

305.5 

334.0 

327.6 

429.4 

438.3 

13 

352.3 

352.1 

356.7 

352.0 

381.1 


454.5 

453.2 

14 

409.1 

410.8 

413.5 

^412.5 

434.7 

436.6 

493.1 

493.3 


a Extensional frequency. 

^Coupled with n = 13. 
c Coupled with n = 9. 

^Maximum frequency with air shaker. 
e Coupled with n = 10. 

^Minimum frequency with electric shaker. 
^Coupled with m = 3, n = 13. 
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TABLE in.- EXPERIMENTAL AND ANALYTICAL INEXTENSIONAL FREQUENCIES 
OF FREE-FREE ELLIPTICAL CYLINDRICAL SHELLS 

(a) £=1.176, e = 0.526 




Frequency, Hz 

, at m = 0 from — 


p 

a Experiment 

Present analysis 

Reference 9 

Reference 10 


Symmetric 

Antisymmetric 

Symmetric 

Antisymmetric 

0 






0 

1 






7.7 x 10-5 

2 

5.6 

5.6 

5.62 

5.68 

5.92 

5.56 

3 

16.1 

16.2 

15.89 

15.89 

16.3 

15.4 

4 

30.9 

30.8 

30.52 

30.52 

31.2 

31.1 

5 

50.1 

50.1 

49.41 

49.41 

50.3 

50.45 

6 

74.8 

74.4 

72.54 

72.54 

73.8 

71.6 

7 

102.4 

102.4 

99.87 

99.87 

101.5 

98.9 

8 

134.6 


131.4 

134.4 

133.5 

133.2 

9 

171.5 

171.7 

167.2 

167.2 

169.8 

168.1 

10 

212.5 

212.8 

207.1 

207.1 

210.3 

204.8 

11 

258.8 

258.4 

251.3 

251.3 

255.2 

251.1 

12 

312.1 


299.6 

299.6 

304.2 

302.8 

13 

363.8 

362.3 

363.0 

352.2 

352.2 

357.5 

351.9 

14 

423.2 


409.0 

409.0 

415.2 

405.8 



Frequency, Hz 

at m = 1 from — 


P 

^Experiment 

Present analysis 

Reference 10 

Reference 11 


Symmetric 

Antisymmetric 

Symmetric 

Antisymmetric 

0 



b 2576.00 


(c) 

0 

1 



b 1886.00 

b 2 148.00 

(c) 

0 

2 

7.7 

7.8 

7.52 

7.55 

6.13 

6.25 

3 

19.1 

19.1 

18.94 

18.94 

16.2 

16.9 

4 

34.6 


34.10 

34.11 

32.0 

33.5 

5 

54.3 

54.3 

53.28 

53.28 

51.5 

50.7 

6 

79.5 

78.9 

76.56 

76.56 

72.6 

73.9 

7 

106.8 

106.9 

104.0 

104.0 

100.0 

101.4 

8 

138.9 


135.6 

135.6 

134.3 

133.0 

9 

174.8 

175.4 

171.4 

171.4 

169.2 

168.9 

10 

215.9 

216.3 

211.4 

211.4 

205.9 

209.0 

11 

261.7 

261.8 

255.6 

255.6 

252.2 

253.3 

12 



304.0 

304.0 

303.9 

301.9 

13 



356.6 

356.6 

353.0 

354.7 

14 

430.1 


413.4 

413.4 

406.9 

411.7 


frequencies for p £ 5 obtained with air shaker. 
b Extensiona.l frequency. 
c Negative eigenvalue. 



TABLE m.- EXPERIMENTAL AND ANALYTICAL INEXTENSIONAL FREQUENCIES 
OF FREE-FREE ELLIPTICAL CYLINDRICAL SHELLS - Continued 


(b) 1.538. e = 0.760 

b 


Frequency, Hz, at m = 0 from — 


Experiment Present analysis 

Symmetric Antisymmetric Symmetric Antisymmetric 

a 2554.00 

0 a 4821.00 

5.1 5.7 5.47 5.89 

15.0 14.9 15.5 15.5 


Reference 9 Reference 10 


1.02 X 10- 4 


98.1 

129.4 

164.9 

d 203.4 

248.8 

295.6 


Experiment 


Frequency, Hz, at m = 1 from - 


Present analysis 


Symmetric Antisymmetric Symmetric Antisymmetric 


Reference 10 Reference 11 


a 2336.00 

a 1694.00 


a 2343.00 


168.2 

f 207.9 

253.1 


a Extensional frequency. 

^Maximum frequency with air shaker. 
c Minimum frequency with electric shaker. 
d Coupled with m = 2, p = 9. 
e Negative eigenvalue. 

^Coupled with m = 2. 



TABLE m.- EXPERIMENTAL AND ANALYTICAL INEXTENSIONAL FREQUENCIES 
OF FREE-FREE ELLIPTICAL CYLINDRICAL SHELLS - Concluded 

(c) — = 2.5, e = 0.916 





Frequency, Hz, at m 

= 0 from — 



p 

Experiment 

Present analysis 

Reference 10 



Symmetric 

Antisymmetric 

Symmetric 

Antisymmetric 


0 



a 2584.00 


0 


1 



0 

a 5346.00 

4.71 x 10- 5 


2 

4.8 


5.13 

6.79 

5.39 


3 

14.0 

13.7 

14.4 

14.4 

16.6 


4 

27.4 

28.5 

28.2 

29.0 

29.3 


5 

b 46.3 

b 45.8 

c 46.9 

46.7 

46.5 

51.5 


6 

69.2 

69.9 

69.2 

69.7 

69.7 


7 

96.4 

96.0 

96.4 

96.2 

103.7 


8 

126.9 

128.3 

127.5 

127.8 

126.7 


9 

163.1 

d 162.7 

163.1 

163.0 

172.9 


10 

202.6 

203.0 

202.9 

203.1 

200.6 


11 

246.5 

247.4 

247.1 

247.0 

258.6 


12 

297.1 

296.8 

295.4 

295.5 

291.8 




Frequency, Hz 

at m = 1 from — 


P 

Experiment 

Present analysis 

Reference 10 

Reference 11 


Symmetric 

Antisymmetric 

Symmetric 

Antisymmetric 

0 



a 3700.00 


(e) 

0 

1 



a 1401.00 

a 2 525.00 

(e) 

0 

2 

7.2 

7.4 

7.11 

7.74 

5.94 

6.76 

3 

b 17.2 

c 18.0 

b 16.2 

17.3 

17.2 

17.5 

18.3 

4 

b 30.9 

b 31.4 

31.7 

32.1 

30.2 

34.4 

5 

b 49.0 

C51.4 

b 49.9 

c 49.5 

50.3 

50.2 

52.6 

55.1 

6 

72.9 

73.1 

73.1 

73.4 

70.7 

80.5 

7 

100.0 

100.4 

100.3 

100.2 

104.8 

110.4 

8 

130.7 

130.8 

131.6 

131.8 

127.7 

144.9 

9 

167.1 

167.3 

167.6 

167.2 

174.0 

184.0 

10 

206.4 

204.8 

207.1 

207.3 

201.6 

227.7 

11 

f 256.6 

g 250.3 

g 260.0 

251.3 

251.2 

259.75 

276.0 

12 

. 

301.4 

300.9 

299.6 

299.7 

292.8 

329.0 


a Extensional frequency. 

^Obtained with air shaker. 
c Obtained with electric shaker. 

^Strongly coupled with m = 1, p = 8. 
e Negative eigenvalue. 

^Coupled with rn = 2. 

^Circumferential nodes (i.e. , for m = 1) at opposite ends of shell. 



TABLE IV.- EXPERIMENTAL AND ANALYTICAL FREQUENCIES OF FREE-FREE ELLIPTICAL CYLINDRICAL SHELLS 



Frequency, Hz, for 

- = 1.176, 
b 

e = 0.526 



Frequency, Hz, for 

2 = 1.538, 
b 

e = 0.760 



Frequency, Hz, for - 

- = 2.5, 
0 

e = 0.916 


Experiment Analysis Experiment Analysis Experiment Analysis Experiment Analysis Experiment Analysis Experiment Analysis Experiment Analysis Experiment Analysis Experiment Analysis 



362.3 
















212.5 

207.1 

363.0 

352.2 

a 508.2 

496.6 

200.7 


332.3 

330.6 


475.6 

148.1 



278.9 


405.4 

212.8 

207.1 

363.8 

352.2 


496.6 

201.7 

206.3 

343.6 

345.3 


476.2 

151.5 

165.6 

260.6 

279.5 


409.6 

215.9 

211.4 


356.6 


500.9 

203.4 

206.3 

344.5 

345.3 

461.5 

476.2 

158.4 

165.6 


279.5 


409.7 

216.3 

211.4 


356.6 


500.9 

207.9 

210.6 

347.9 

345.4 


477.7 

163.2 

166.1 

296.8 

295.4 


411.7 

241.4 

244.1 

375.9 

373.8 



208.4 

210.6 


345.4 


477.7 


166.1 

297.1 

295.5 


412.3 

241.8 

244.1 

376.6 

373.8 



211.9 

213.9 

349.5 

351.7 


484.5 

167.1 

167.2 

300.9 

299.6 

423.0 

431.9 

250.6 

246.5 

390.6 

381.5 



218.4 

213.9 


351.7 


484.5 

167.3 

167.6 

301.4 

299.7 


432.2 


246.5 

391.0 

381.5 




214.0 


356.0 



202.6 

202.9 


327.9 


476.7 

258.4 

251.3 


400.6 




214.0 


356.0 



203.0 

203.1 


327.9 


476.9 

258.8 

251.3 


400.6 



248.8 

250.5 

362.2 

366.6 



204.8 

207.1 


328.4 


478.7 

261.7 

255.6 


401.6 



248.9 

250.5 

367.2 

366.6 



206.4 

207.3 


328.4 


479.0 

261.8 

255.6 


401.6 



253.1 

254.8 

380.8 

383.8 




217.6 


340.6 



a 279.7 

274.4 

423.1 

409.0 



253.3 

254.8 


383.8 



209.1 

217.6 


340.6 




274.4 

423.2 

409.0 



259.2 

264.3 

399.4 

410.2 




219.2 


347.8 



304.5 

292.1 

430.1 

413.4 



265.2 

264.3 

401.0 

410.2 



246.5 

219.2 


347.9 




292.1 


413.4 




268.8 

403.9 

410.5 



247.4 

247.0 


352.2 



312.1 

299.6 

445.8 

434.8 




268.8 


410.5 



250.3 

247.1 


352.3 




299.6 

450.9 

434.8 



295.6 

298.9 

418.0 

411.9 



255.1 

251.2 

360.6 

356.5 




304.0 


444.6 



296.3 

298.9 


411.9 



256.6 

251.3 

363.0 

356.5 



295.7 

304.0 


444.6 



301.4 

303.3 


414.7 




265.3 


398.6 



297.9 

304.4 , 

465.0 

452.8 



301.7 

303.3 


414.7 




265.3 


398.7 



306.7 

304.4 

466.3 

452.8 



308.5 

311.5 

440.4 

438.1 



260.0 

265.4 


399.4 



340.3 

334.8 

476.1 

484.2 




311.5 


438.1 




265.4 


399.4 



| 341.4 

j 






• 











343.7 

334.8 

a 491.0 

484.2 


1 I 

323.4 

330.6 


475.6 

J 

1 

260.4 

278.9 


405.3 




a Obscure circumferential nodal pattern. 


CO 


TABLE V.- CONVERGENCE CHARACTERISTICS FOR 
VARIOUS CIRCUMFERENTIAL MODAL 
APPROXIMATIONS WITH A SINGLE 
LONGITUDINAL TERM 



e = 0.916; 


symmetric 


p 

Number of circumferential terms 

30 

32 

34 

35 



m = 0 



0 





2 

5.126 

5.126 

5.126 

5.126 

4 

28.23 

28.23 

28.23 

28.23 

6 

69.19 

69.19 

69.19 

69.19 

8 

127.5 

127.5 

127.5 

127.5 

10 

203.0 

202.9 

202.9 

202.9 

12 

295.5 

295.4 

295.3 

295.3 

'U 

407.2 

405.4 

404.8 

404.7 



m = 1 



0 

a 0 

a 0 

a 0 

a 0 

2 

7.106 

7.106 

7.106 

7.106 

4 

31.70 

31.70 

31.70 

31.70 

6 

73.07 

73.07 

73.07 

73.07 

8 

131.6 

131.6 

131.6 

131.6 

10 

207.1 

207.1 

207.1 

207.1 

12 

299.8 

299.6 

299.6 

299.6 

14 

411.4 

409.6 

409.1 

409.0 


a u-mode. 



TABLE V.- CONVERGENCE CHARACTERISTICS FOR 
VARIOUS MODAL APPROXIMATIONS - Concluded 


= 2.5; e = 0.916; symmetric 


3519.0 

908.8 

5209.0 

7034.0 


Number of circumferential terms 
31 32 j 33 34 

m = 2 


3519.0 

908.5 

5209.0 

7034.0 


3519.0 

908.0 

5209.0 

7034.0 


3519.0 

906.8 

5209.0 

7034.0 


3519.0 

5209.0 

7034.0 


3519.0 

911.5 

5209.0 

7034.0 


167.6 

536.2 

221.3 

282.7 

343.5 

434.8 

581.6 


167.6 

534.9 

221.3 

282.7 

343.3 

433.7 

577.1 


221.3 

282.7 

343.3 

533.9 

433.0 

574.0 


167.6 

221.3 

282.7 

343.2 

533.2 

432.7 

571.9 


167.6 

221.3 

282.7 

343.2 

532.7 

432.5 

570.6 

902.2 


167.6 

221.3 

282.7 

343.2 

532.4 

432.4 

569.7 


5415.0 


1400.0 

861.8 

269.9 

333.8 

405.6 

485.1 

664.2 

584.0 


1396.0 

1396.0 

1395.0 

1394.0 

1394.0 

5415.0 

5415.0 

5415.0 

5415.0 

5415.0 

8521.0 

8521.0 

8521.0 

851.5 

8521.0 




8521.0 



858.3 

269.9 

333.7 

405.4 

483.9 

660.3 

580.6 


855.6 

269.9 

333.7 

405.2 

483.2 

657.4 

578.3 


853.4 

269.9 

333.7 

405.1 

482.8 

655.2 

576.6 


269.9 

333.7 

405.1 

482.6 

653.6 

575.6 



TABLE VI.- RIGID-BODY (ZERO FREQUENCY) MODES 










m 

P 

Mode 

e = 0 

e = 0.526 

e = 0.760 

e = 0.916 


0 

Symmetric 


— 

— 

— 

0 


Antisymmetric 

V 

V 

V 

V 






‘ ■ 


1 

Symmetric 

v or w 

V 

w 

w 



Antisymmetric 

v or w 

V 

V 

V 


0 

Symmetric 

u 

u 

u 

u 

1 


Antisymmetric 


— 

— 

— 





■ 

— 


1 

Symmetric 

u 

V 

V 

a v 



Antisymmetric 

u 

u 

u 

u 

Number of 

J~Longitudinal 

1 or 4 

4 

4 

2 

terms 

^Circumferential 

1 

13 

16 

32 


a Actual frequency, 1.01 Hz. 


34 



TABLE VII.- ANALYTICAL FREQUENCIES FOR FREELY SUPPORTED 

CIRCULAR SHELL (MODEL 1) 


n 

0 

1 

2 

3 


m = 1 

2537.0 

1565.0 
894.1 
529.8 


Frequency, Hz 

m = 2 

2594.0 

2309.0 

1782.0 

1315.0 



TABLE Vni.- ANALYTICAL FREQUENCIES FOR FREELY SUPPORTED 

ELLIPTICAL SHELLS 





Frequency, Hz 



p 

e = 0.526 (model 2) 

j e = 0.760 (model 3) 

e = 0.916 (model 4) 


Symmetric 

| Antisymmetric 

| Symmetric 

j Antisymmetric 

| Symmetric 

| Antisymmetric 




m 

= i 



0 

2550.0 


2612.0 


2723.0 


1 

1440.0 

1686.0 

1238.0 

1856.0 

931.1 

2039.0 

2 

876.6 

888.9 

785.2 

858.5 

555.6 

748.1 

3 

524.1 

524.2 

491.1 

492.4 



4 

335.5 

335.5 

319.4 

319.4 



5 

234.3 

234.3 


226.9 


392.4 

6 

184.2 

184.2 



106.8 

106.8 

258.2 

7 

157.1 

157.0 

138.5 

138.5 

107.0 

388.4 

107.7 

151.6 

8 

160.2 

160.2 

140.1 

178.3 

140.1 

178.3 

148.2 

257.0 

148.2 

9 

189.8 

189.8 

184.1 

226.9 

184.1 

151.7 

194.8 

194.9 

10 

221.9 

221.9 

223.9 

223.9 

206.5 

206.6 

11 

261.9 

261.9 

263.6 

263.6 

259.3 

259.4 

12 

308.1 

308.1 

307.3 

307.3 

310.2 

310.5 

13 

359.5 

359.5 

359.4 

359.4 

351.4 

351.9 

14 

415.6 

415.6 

417.1 

417.1 

413.3 

413.5 




m 

= 2 



0 

2726.0 


3248.0 


4155.0 


1 

208 3.0 

2607.0 

18 76.0 

3188.0 

1694.0 

4081.0 

2 

1696.0 

1782.0 

1419.0 

1775.0 

1035.0 


3 

1291.0 

1297.0 

1145.0 

1214.0 

793.2 

989.0 

4 

956.5 

956.7 

886.1 

893.0 


1675.0 

5 

718.8 

718.8 

677.0 

677.4 



6 

555.6 

555.6 

531.3 

531.6 



7 

446.6 

446.6 

441.5 




8 

381.8 

381.8 



659.7 


9 



268.3 

268.3 

202.4 

202.4 

10 

310.6 

310.6 

268.1 

328.2 

268.1 

202.5 

258.5 

202.5 

258.5 

11 

312.4 

347.5 

312.4 

347.5 

325.4 

325.4 

259.2 

259.2 

12 

359.4 

359.4 

382.3 

328.2 

382.3 

323.8 

398.0 

323.8 

398.2 

700.3 

13 

401.5 

401.5 

393.4 

393.4 

441.5 

393.2 

550.4 

393.4 

14 

448.4 

448.4 

454.7 

454.7 

466.3 

467.2 
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Figure 1.- View of free-free elliptical cylindrical shell and circumferential track. 




Figure 2.- Schematic of simplified test apparatus. 
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Figure 3.- Closeup of movable inductance probe assembly 



Longitudinal 


p * 3 


p ■ 4 


Symmetric Anti-symmetric 

Circumferential 

Figure 4.- Mode-shape classifications for normal (or radial) displacements w of 
free-free elliptical cylindrical shell. For the circular shell, p = n where n 
is the number of circumferential waves. 
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Figure 5.- Continued 
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Figure 5.- Continued. 
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Figure 5.- Continued. 
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Figure 7.- Circumferential mode-shape correlation for free-free circular 
cylindrical shell (model 1; e = 0). 
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Analysis 


Experiment 



0 2 A .6 .8 1-0 

s/s o 

(a) Inextensional symmetric mode, m = 0. 

Figure 9.- Circumferential mode-shape correlation for model 3 (e = 0.760) 
Experimental mode shapes taken from figure 5. 
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(f) Antisymmetric modes involving m = 2. Frequencies about 210 Hz. 

Figure 9.- Continued. 
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Experiment 
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(g) Symmetric modes involving m = 2. Frequencies about 250 Hz. 

Figure 9.- Continued. 








(k) Symmetric modes involving m = 3. Frequencies about 350 Hz 

Figure 9.- Continued. 
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Analysis 


Experiment 
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(n) Antisymmetric modes involving m = 3. Frequencies about 400 Hz. 

Figure 9.- Continued. 
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Normalized w-displacement 


Analysis 


Experiment 


f 3 | | - 411.9 Hz 



(p) Antisymmetric modes involving m = 3. Frequencies about 425 Hz. 

Figure 9.- Continued. 
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Experiment 


f 2,5 * 475.6 Hz 




(r) Antisymmetric modes involving m = 3. Frequencies about 475 Hz. 

Figure 9.- Concluded. 
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Figure 11.- Longitudinal mode-shape correlation for model 4 (e = 0.916) 
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Figure 12.- Effect of various circumferential modal approximations on mode shapes 

of model 4 (e = 0.916). 
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Figure 15.- Variation of frequencies with mode numbers for freely supported shells 
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